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ABSTRACT 


A  systematic  theoretical  development  is  presented  for  sound 
absorption  in  dilute  chemically  reacting  systems  when  heat  effects 
and  other  dissipative  mechanisms,  e.g.  diffusion  and  viscous  damping, 
are  ignored. 

A  derivation  is  presented  for  the  isothermal  frequency  depen¬ 
dent  compressibility  for  such  a  system  and  the  relation  of  this 
quantity  to  sound  attenuation  and  the  sound  absorption  spectrum. 

The  method  is  applied  to  several  illustrative  reaction 
mechanisms  of  interest.  Particular  attention  is  placed  on  the  sim¬ 
plifications  that  arise  when  there  is  a  single  slow  step  in  the 
reaction  mechanism.  The  general  case  of  coupled  chemical  reactions 
is  treated  and  a  prescription  presented  for  detei-mining  the  sound 
absorption  coefficient. 

The  final  section  is  concerned  with  sound  absorption  from 
chemical  systems  in  steady  states  far  from  chemical  equilibrium. 

For  such  systems  the  possibility  exists  that  composition  fluctuations 
will  behave  as  chemical  waves  which  will  be  exhibited  by  a  distinctiv 
splitting  in  the  sound  absorption  spectrum,  i 


I .  Introduction 

Sound  absorption  in  chemically  reactive  fluids  is  a  subject  of 
interest  to  workers  in  a  wide  variety  of  fields,  including  chemistry, 
biology,  and  fluid  mechanics.  The  subject  also  has  some  practical 
importance.  In  particular,  sound  absorption  at  relatively  low  fre¬ 
quency  (below  a  few  hundred  kilohertz)  in  the  ocean  is  markedly 
influenced  by  chemical  relaxation  between  the  prodiguous  number  of 
different  species  (mostly  ionic)  that  are  present  in  sea  water. 
Understandably,  the  subject  has  attracted  the  attention  of  the  Navy 
and  the  Defense  Advanced  Research  Projects  Agency  (DARPA;  who  have 
supported  work  in  this  area  over  the  years.  Substantial  progress  has 
been  made  in  understanding  the  rather  complicated  chemistry  and 

resulting  sound  absorption  in  sea  water.  However,  some  problems 

f  1 ' 

remain,  notably  the  absorption  at  Ik-Hz  attributable^  '  to  boric 
acid  in  sea  water. 

Consideration  of  this  oroblem  and  the  methods  of  analysis  em¬ 
ployed  to  interpret  the  sound  propagation  measurements  indicates  a 
need  for  a  systematic  development  of  ^-heory  of  sound  absorption 

in  chemically  reacting  systems.  While  theory  is  available  in  a 
(2-5) 

number  of  sources  ,  the  existing  developments  are  either  overly 
formal  or  difficult  to  penetrate.  Moreover,  the  existing  developments 
do  not  describe  the  approximations  needed  to  reduce  tf',e  general  case 
of  the  influence  of  coupled  reactions  in  solution  on  sound  attenuation 
to  the  frequently  eni.ountcred  -"^se  where  a  sinnle  step  in  the  reaction 
mechanism  is  slow  and  hence  rate  determining. 
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The  purpose  of  this  technical  note  is  to  provide  a  clear  and 
direct  theoretical  development  for  sound  absorption  arising  from 
chemical  relaxation.  Hopefully  this  development  will  be  useful,  not 
only  for  those  working  in  sound  attenuation  in  the  oceans  but  also 
for  others  working  on  sound  absorption  in  chemistry  and  biology. 

The  treatment  presented  is  restricted  to  the  circumstances 
most  frequently  encountered.  It  is  assumed  that  dissipative  trans¬ 
port  processes  other  than  chemical  reaction,  e.g.  diffusion,  thermal 
conduction,  and  viscous  damping,  can  be  ignored,  which  is  usually  the 
case  at  low  frequencies.  Thermal  (enthalapic)  effects  due  to  heats 
of  reaction  are  also  ignored  so  the  treatment  is  for  isothermal 
conditions.  Moreover,  the  analysis  is  restricted  to  dilute  solutions 
so  that  all  activity  or  concentration  effects  are  ignored.  All  of 
these  restrictive  assumptions  except  for  the  last  one  (low  concentra- 
tTon)  can  easily  be  relaxed. 

In  Section  II,  the  pertinent  solution  thermodynamics  is  pre¬ 
sented.  Section  III  presents  the  oertinent  relaxation  equations  for 
this  problem,  and  in  Section  IV  the  determined  frequency  dependent 
isothermal  compressibility  is  related  to  sound  absorption.  Section  V 
is  devoted  to  demonstrating  for  several  mechanistic  examples  of 
interest  how  the  method  may  be  employed  to  determine  the  sound 
absorption  spectrum.  Section  VI  is  concerned  with  a  description  of 
the  general  case  that  may  be  encountered  in  coupled  chemical  reactions 
where  the  assumption  of  a  single  slow  step  is  not  valid.  The  analysis 
this  section  would  be  required  if  one  wished  to  model  a  complex 
cnemically  reacting  system,  e.g.  sea  water,  comnrehensively  and  not 
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divide  the  problem  so  to  sneak,  piecewise,  where  each  frequency 
range  is  associated  with  a  single  dominant  relaxation  step.  Tne 
final  section  is  concerned  with  sound  absorption  from  chemically 
reacting  systems  in  a  steady  state,  far  from  equilibrium,  where 
the  fluctuations  in  the  chemical  reaction  behave  as  "i,liemical  waves" 
and  influence  the  attenuation  of  the  sound  waves  accordingly.  This 
last  section  contains  material  which  should  be  of  interest  to  the 
many  workers  who  are  currently  devoting  attention  to  non-linear 
chemically  reacting  systems. 
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II.  Solution  Thermodynamics 

For  0  solution  containing  several  components,  the  Gibbs  free  energy 

^  1*  c* 

o  1  S 

(2  1 ) 

'S  G-  ^  ^  4V  4  ,5  ^/v\c 


where  S  in  the  entropy,  T  the  absolute  temperature,  V  the  volume,  and 
the  pressure.  Here  iii  is  the  chemical  potential  of  solute  species  i 
and^'^i  the  number  of  moles  of  species  i;  and/^Q  denote  resoectively 
the  chemical  potential  and  moles  of  solvent,  assumed  inert. 

We  assume  that  the  solute  species  are  undergoing  r  chemical 
reactions  of  the  form 

X  ^2.2) 

...,r 


•f"  h 

with  the  convention  that  the  stochiometric  coefficients  in  the 
reaction,  ,  are  positive  for  products  and  negative  for  reactants. 
We  define  the  affinity  of  the  reaction  Aa  as 


A<^ 


(2.3) 
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•f  K 

and  the  progress  variable  for  the  reaction  <sxa  as 

(oc)  ^  (2.4) 

fO*  -  9;  . 


Since  the  interconversion  of  species  can  only  occur  by  chemical  reaction 


z: 


(2.51 


and  Eq  (2.1)  may  be  expressed  as 


(2.6) 


The  solvent  is  assumed  not  to  participate  in  any  of  the  reactions. 

At  chemical  equilibrium  (cGl-y,^  =  o  which  yields  the  condition 

(2.7) 

Ao<--A.  CT,t>,  A')  =0 

which  determines  the  equilibrium  values  of  the  variables  a°  =  (a°^  .  . 
•°^).  These  values  will  deoend  upon  T  and  ,  so  =  \^^(T,p). 

We  shall  assume  that  the  species  in  solution  that  undergo  chemical 
reaction  exhibit  ideal  behavior  so  that 
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1 


I 


(2.9) 


where  x^-  :  n.j/no  is  the  mole  fraction  of  species  i.  It  follows  that 
we  are  restricting  consideration  to  dilute  solutions. 

We  shall  be  concerned  with  isothermal  relaxation  processes,  i.e. 
temperature  and  enthalpy  effects  will  be  neglected.  The  formulation 
can  easily  be  generalized  to  include  these  processes. 

The  differential  change  in  the  voume  V  is  given  by 

(2.10) 


The  equilibrium  isothermal  compressibility  <-|-(0)  is  defined  by 


(2.11) 


and  the  "frozen"  isothermal  comoressibil  ity  •Cy('^)  is  defined  by 


v  w  y 


(2.12) 


8 


III  Relaxation  Behavior 

We  imagine  a  pressure  wave  of  the  form 

(3.1) 

imoosed  on  the  reacting  mixture.  This  disturbance  will  instantaneously 
perturb  the  chemical  environment.  The  soecies  will  interconvert 
in  an  effort  to  attain  the  new,  momentarily  established,  local 
equilibrium  condition  A^CT,  p+op,  a+6>0  =  0.  For  small  disturbances 
the  changed  equilibrium  is  gi'-en  oy 


The  chemical  species  will  not  be  able  to  follow  the  pressure 
disturbance  instantaneously,  and  this  v/ill  cause  relaxation.  The 
relaxation  behavior  is  determined  by  the  kinetic  equations  which  are 
governed  by  the  linear  laws 
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(3.3) 


The  factor  V~^  included  on  the  left  hand  side  takes  into 
account  that  the  progress  variables  as  defined  here  are  extensive 

I 

proDortios;  the  kinetic  coefficirnt-  >  p  are  intensive  quantities. 
This  expansion  may  be  written  in  more  compact  vector  notation 


i 

V 


(3.4) 


13.5) 


(3.6) 


In  vector  notation  Eq.  (2.10)  is 
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(3.7) 


where  use  has  been  made  of  the  Maxwell  relation  which  follows  from 
Eq.  (2.6) 

(3.8) 


Jt 


iXV 


— ■*  I  -  lT 


At  equilibrium  =  0  and  Eq.  (3.4)  yields  the  relation 


(3.9) 


Thus,  the  compressibil ity  relation,  Eq.  (2.13),  can  as  well  be  expressed 
as 

(3.10) 


X-r(o)  = 


a  relation  that  will  prove  useful  shortly. 

The  imposed  pressure  will  induce  a  steady  state  response  in  the 
volume  and  progress  variables  given  by 

^  A  ^  (3.11) 


n 


respectively.  The  steady  state  resoonse  for  isothermal  channes  in 
volume  is 


(3.12) 


A 

We  may  solve  for^^  by  substitution  Eq.  (3.11)  i.i  Eq.  (3.12);  the 
result  is 


=  _v  I  xiol  +  L-^l]  .  L- J;  H  . 


(3.13) 


Before  proceeding,  it  is  useful  to  write  Eq.  (3,13)  in  terms  of 
more  familiar  rate  expressions  by  defining  the  kinetic  matrix 


M 


(3.14) 


then  Eq.  (3.13)  takes  the  form 


_\ 


(3.15) 
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In  terms  of  time  Eq.  (3.15)  corresponds  to 


,  -h.f)  -  VH-Sc 


(3.16) 


t 


—  ^ 


where  5£  are  "imposed"  concentration  differences  induced  by  the  pressure 
wave  perturbation 


(3.17) 


We  continue  by  eliminating  5X  between  Eqs.  (3.13)  and  (3.12) 


thus  obtaining  a  relation  between  5V  and 


(3.18) 


^  /\ 

The  quantity  'iV/'^^can  be  related  to  the  frequency  dependent 
isothermal  compressibility 


(3.19) 
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Then  Eq.  (3.18)  may  be  written  as 


(3.20) 


X-j-  ~  ') 


+  !C 


A1  ternatively,  Eq.  (3.10)  may  be  employed  to  express ^j(‘J)  as 

(3.21) 

YjL^)  ~  +  , 


In  the  limit  '^^0,  <j(a))^<.|.(0) ,  and  in  the  limit  ,  <j(a) )-^<j(”) . 

This  expression  for  <y(<i))  is  the  central  starting  point  for  the 
investigation  of  sound  absorption  in  chemically  reacting  fluids.  It  is 
important  to  bear  in  mind  the  restrictive  assumptions  under  which  this 
result  has  been  obtained.  In  particular,  temperature  effects  have  been 
ignored  and  local  equilibrium  conditions  have  been  assumed.  In 
Section  VI  we  will  indicate  how  this  expression  can  most  easily  be 
evaluated  for  complex  chemically  reacting  systems  of  interest. 
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IV  Relation  to  Sound  Attenuation 

The  isothermal  sound  speed  Cj  is  given  by 


(4.1) 


and  the  pressure  satisfies  the  wave  equation 


a 


(4.2) 


when  viscous,  diffusive,  and  heat  dissipation  effects  are  ignored. 
For  a  pressure  wave  of  the  form  Eq.  (3.1)  one  obtains  the  dispersion 
relation 


Cd 


(4.3) 


For  the  case  of  chemical  relaxation  the  corresponding  dispersion 
relation  is 


(4.4) 


and  the  frequency  dependent  compressibility  leads  to  both  sound 
absorption  and  velocity  dispersion.  For  the  case  of  an  impressed 
pressure  wave  disturbance  of  real  frequency  w,  the  complex  frequency 
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dependency  in  Cj('^)  will  force  R  to  become  complex. 
Let 

k  =  L'-k'' 


(4.5) 


and  seek  a  solution  to  Eq.  (4.4)  for  the  imaginary  part, 
express  <y(tLi)  in  real  and  imaginary  carts 


-.Wylc^-)  -i  V'-rto) 


one  finds 


If  we 


(4.6) 


(4.7) 


Generally,  one  is  justified  in  assuming  that  the  characteristic  length 

-1  i'  -1 

for  attenuation  (|^")  is  large  compared  to  the  wave  length  ( |^  )  so 

p 

that  terms  of  order  (V/V)  in  Eq.  (4.7)  may  be  ignored.  Under  these 
circumstances  one  finds 


(4.8) 


and 
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(4.9) 


The  speed  of  sound 


is  simply  related  to  according  to 


(4.10) 


The  resulting  expression  for  the  pressure  wave 


(4.11) 


clearly  identifies  as  the  attenuation  in  the  medium, 
per  wave  length  a(a))  is  defined 

o{(U)  -  ^  I ^  , 

Consider  the  simple  case  of  single  relaxation  time  where 


'X^fu-y-.y^roc',  +  A 


The  adsorption 


1 


(4.12) 


(4.13) 


then 
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^  (t^)  4-  — ^ 


(4.14) 


and 


/  4 


(4.15) 


so  that 


I  +  (U7:f 


(4.16) 


and 


gj  2: 

/  f  fioz)^ 


(4.17) 


In  the  general  case  of  many  relaxation  times 

)4-  +  2r  ^ 


(4.18) 
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the  calculation  of  ci((ij)  and  Cj{(,})  presents  an  algebraic  task  that,  in 
general,  will  require  numerical  calculation.  However,  for  many  cases 
of  interest  the  dispersion  in  velocity  [Ag^/<y(«>)]  will  be  quite  small 
and  a(co)  can  be  approximated  by 

0(10)  i  Z  ^  , 

>7/0^)  ^  l-h(cJli)  > 


(4.19) 
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V  Some  Examples 

First  consider  the  case  where  there  is  no  coupling  between  the 
reactions.  Then 


L.  “  ^  j  t/1,  = 


(5.1) 


and  it  follows  from  Eq.  (3.20)  that 


(5.2) 


rpAvN 

<?>•«  >'T;P'  i  I  +  t  ^ 


Note  that  this  expression  has  the  form  of  Eq.  (4.18). 

For  a  single  reaction,  the  result  of  Eq.  (3.20)  or  (5.2)  is 

(5.3) 


^  (u')- 


I 

I  +A{^Z 


V . A .  First  Order  Equilibria 

Consider  a  simple  unimolecular  eouilibrium 

k,  (5.4) 
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For  this  mechanism 

(5,5) 


or  since  the  total  number  of  moles  is  constant 


(5.6) 


where/ri^  andA\0  are  the  number  of  moles  in  the  volume  V.  The  progress 
variable  is  simply 


(5.7) 


The  derivatives  required  for  the  evaluation  of  Eq.  (5.3)  are  easily 
found. 

(5.8) 


3^ 

B'X  Jt 


l^\ 


and 


^  RT 


+ 


A)V 


m 


tj 


(5.9) 


21 


with/»^  =  rr?^  total  number  of  moles  in  volume  V  in  the 

solution. 

It  should  be  noted  that 


and  since  the  equilibrium  constant  K  is 

K  =  -eitp  [- 


we  may  find  AvrfromjjC  according  to 


^vr  =  -  n 

V  »\> 


The  relaxation  time  must  be  determined  from  the  kinetic  equations 
which  we  write  in  terms  of  the  moles  of  species  in  the  volume  V  instead 
of  the  concentrations.  This  is  required  because  the  progress  variables 
have  been  defined  as  extensive  quantities. 

-  -  le, 

^  -  I?-' 
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Note  that  at  equilibrium 

(5.11) 

k,  »’a  -  k-."!  ov  Ji?  -  =  K, 

K  k-' 

SO  that  in  terms  of  deviations  from  equilibrium  Eq.  (5.10)  is 

(5.12) 

d  ^  i?-! 

fc 

jjm,  ^  -  fc.,  r>'6 , 

In  terms  of  the  progress  variable,  Eq.  (5.7)  one  obtains  a  single 
equation 

(5.13) 

^  +  t-.) 

t 


■ 
I 


which  immediately  identifies  the  relaxation  time  as'^"^  =  (j^^ 
In  sum^  <rj(v}  is 


y^oo)  + 


[- x;x;  1  ' 

(Xa+'< 


(5.14) 


whe^'e  ^  =  V//K\. 
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V . B .  Coupled  First  Order  Equilibrium  With  a  Slow  Step 
We  next  consider  the  coupled  reaction  scheme 


c 


(5.19) 


which  has  the  following  kinetic  equations  (in  terms  of  moles  in  the 
volume  V) 

=  -  If,  (V\;^  -i-  Is-, 


at 


and  two  progress  variables 


(5.18) 


In  terms  of  the  progress  variables  one  has 


At 

(5.19) 

d 

--  n.  -  . 

d  t 

J 
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We  consider  the  special  case  where  t?:c  second  reaction  is  slow  compared 
to  the  first.  In  this  case  one  may  set  =  0  in  Eq.  (5.19)  and 
obtain  the  relation 


(5.20) 


and  hence 


(5.21) 


where  we  have  introduced  the  equilibrium  constants 


K 


1 


"'a 


(5.22) 


Thus  the  slow  rela.xation  is  characterized  by  the  relaxation  time 


)  -VK,  i 


(5.23) 
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The  effect  of  assuming  the  rapid  equilibration  of  the  first  step 
[and  the  steady  state  approximation  (JSJi/Ji)]  is  that  the  coupled 
mechanism  may  be  effectively  reduced  to  a  single  step  reaction  mechanism 

(5.24) 


A 


c 


where 


I  VK, 


with  kinetics 

4  k-t 


(5.25) 


(5.26) 


or,  in  terms  of  progress  variables 

41^  r  -  (  ^  ^ 

<?((r 


(5.27) 


which  is  identical  to  Eq.  (5.21). 

We  note  that  for  the  effective  single  step  process,  Eq.  (5.24) 
the  equilibrium  constant  is 


J 
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K.K, 

1  +  K)  • 


(5.28) 


We  can  determine  the  effective  volume  change  for  the  reaction  in 
the  normal  way 


(5.29) 


which  immediately  yields 


1 


(5.30) 


where 


^  H 


iTa  -  U? 


We  may  now  express  the  frequency  dependent  isothermal 
compressibility  for  this  effective  one-step  process  as 

(A 1  - _ ! — 

ITu-  L  x^V+xJ  i  I  4  xw-x 
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in  analogy  to  Eq.  (5.14)^  where  x  is  given  by  Eq.  (5,23)  and  av*  by 
Eq.  (5.32).  It  follows  from  Eq.  (4.12)  that  in  the  limit  of  low 
dispersion  the  absorption  coefficient  is 


TT  (  ALT^) 


L  x;,  +  x: 


cj  '2: 


\  j  4 


(5.31) 


which  may  be  written  in  the  alternative  form 


Rr 


t 

C 


Cf  i  tl 


(5.32) 

1  V 


where  etc.  are  the  equilibrium  concentration.  Using  the 

definition  of  x  Eq.  (5.23)  and,  according  to  Eq.  (5.28) 


this  expression  may  be  simplifipd  to 


(p<  (u)  - 


li  ^  tl  ur' 


(5.33) 
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For  purposes  of  later  comparison  to  the  work  of  others  we  note  that 

(5.34) 

c: 

I  -f 


so  that 


(5.35) 


F  ^  I  +  Ki 


)  4-  (Wt)’- 


V.C.  More  Complicated  Coupled  Equilibria  with  a  Slow  Step 

We  note  that  more  complicated  reaction  schemes  that  contain  a 
slow  step  may  be  treated  by  the  same  procedure  as  described  in  the 
previous  section.  As  a  final  example  we  consider  the  scheme  extensively 
discussed  by  Eigen  and  de  Maeyer^^^  and  employed  by  Mellen  et  al_.  in 
their  investigations  of  the  low  frequency  sound  absorption  in  sea 
viater^‘\  The  scheme  is 

le,  V.  (5-36) 

A  +  ^  7~zf 

fe-, 


The  kinetic  equations  are 

(5.37) 

dA-  .  _v,A(2,  +  W>C.  =  i-? 


^  4,  Ali  -  (k-.fWU 

i£  -  4  k  r  - 

ll  -  ^ 


where  we  have  simplified  notation  by  writing  ^  etc.  In  terms  of 
small  deviations  from  equilibrium,  we  have 


• -k(A«S44?>.5M 

+k-» 

‘  ‘ill 

Ji  b 

-b  k-»-  ^  ^ 

+  V,  Ic 

-  k.. 

For  this  reaction  the  progress  variables  are 

(5.39 

— 5P-.  fc-v-’(5>,-iA0 

where  the  factor  V  is  included  in  order  to  relate  the  change  in  con¬ 
centrations  in  the  volume  V  to  the  change  in  the  extensive  progress 
variables.  So  we  obtain 


+  ki  S  ^ 


(5.40) 
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which  is  precisely  in  the  same  form  as  Eq.  (5.19),  if  we  make  the 
identification 

(5.41) 

-  k,  Ao  . 


If  again  we  assume  the  second  step  is  rate  limiting  and  set  (e)^^i/cit)  = 

we  find,  exactly  as  in  the  previous  example,  that 


-f'  =  V-i  [  1  + 

where 

A 

K,  *  =  K,  ( A.  +1.)  , 


(5.42) 


(5.43) 


0, 


We  may  again  view  the  rapid  first  step  as  establishing  an  equilibrium 
that  leads  to  an  effective  single  step  mechanism  of  the  form 


with,  now,  [c.f.  Eq.  (5.25)] 


(5.45 
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and  associated  equilibrium  constant  [c.f.  Eq.  (5.28)] 

^  (5.46) 

It  follows  that  for  this  effective  one-step  reaction  the  change  in 
volume  is 

A  (5.47) 

t\\r*  AOi  +  /0+  It,). 

The  remainder  of  the  analysis  for  the  absorption  spectrum  follows  the 
preceeding  example  exactly.  We  obtain  the  result 


KT  I  f  I  +  ioT 


with  T  given  by  Eq.  (5.42),  This  result  is  identical  to  that  obtained 
(!) 

by  Mel  1  en  ^ 

It  has  been  correctly  pointed  out  that  this  method  may  be  applied  to 
other  coupled  reaction  mechanism  (provided,  presumably,  that  unstable 
states  are  not  present)  where  a  single  slow  step  is  rate  limiting. 

This  result  is  presented  next. 
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Consider  the  general  case  where  the  slow  step  is  first  order 
between  species  Y  and  Z.  Each  of  these  species  may  be  involved  either 
directly  or  indirectly,  through  coupling  to  other  species,  in  reactions 
that  are  rapid  compared  to  the  interconversion, 

k,  (5.48) 

^  2-  . 


If  one  writes  the  entire  reaction  mechanisims  in  terms  of  the  pertinent 
rate  equations  for  the  progress  variable  and  sets  a  11  =  0 

except  for  the  variable  S  characterizing  the  slow  step  Eq.  (5.43),  one 
arrives  at  the  effective  rate  equation 


/ 


(5.49) 


which  corresponds  to  the  reaction  scheme: 


(6.50) 


The  rates  and  will  be  found  as  a  function  of  the  rate  co¬ 
efficients,  equilibrium  constants,  and  equilibrium  concentrations  of 
other  species  in  the  mechanisim.  The  change  in  volume  AV*  is  given  by 


_L  p 

^  y'v 


(5.51) 
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where 


K 


.af/t-!] 


(5.52) 


and  the  absorption  spectrum  is  simply  expressed  as 


Ri\r  L 


(5.53) 


14  t.,]  ^ 


Thus  we  have  presented  an  explicit  method  for  determining  the 
absorption  coefficient  for  the  case  of  a  single  slow  first  order  step 
in  a  complex  coupled  chemical  reacting  system.  The  method  is  easily 
generalized  to  the  case  where  the  slow  step  involves  the  species  Y  and 
Z  in  higher  order  reactions. 
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VI  Analysis  of  the  General  Case  of  Coupled  Reactions 

We  return  to  Eq.  (3.20),  the  general  expression  for  kj(u))  for  a 
coupled  chemical  system.  We  shall  obtain  an  expression  for  <j(a))  and 
hency  a(aj)  in  terms  of  the  eigenvectors  and  eigenvalues  x  of  the 
kinetic  matrix  M.  In  general  M  is  not  a  symmetric  matrix  (reflecting 
the  fact  that  equilibrium  constants  are  not  unity)  so  it  is  necessary 
to  be  concerned  with  both  right  and  left  handed  eigenvalues  of  M: 

M  ’  A.  t. 


and 


(6.2) 


al  ternatively 


(6.3) 


These  eigenvectors  are  constructed  to  be  orthonormal 


(6.4) 
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and  form  a  complete  set  in  the  vector  space  of  the  r  chemical  reactions. 
Thus  we  may  expand  any  vector  t  as 

(6.5) 

f  =  Z  a.  +>. 

(Vv 


and  the  coefficient  ^J^will  be  given  by 

(6.6) 

fl  /vi  -  (  tv,  )  *  i  • 


Alternatively  the  vector  ^  may  be  expanded  as 


unT 


f  -  I  Up 


'icient  |j 


and  the  coefficient  U  will  be  given  by 


w  •  ly- . 


(6.7) 


(6.8) 


We  shall  expand  the  vectors  (r^  .j7)  according  to  Eq.  (6.5) 
and  the  vectorvT  according  to  Eq.  (6.7)  and  substitute  these  expressions 
into  Eq.  (3.20).  If  we  make  use  of  Eqs.  (6.1)  and  (6.4)  we  obtain 
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where 


It  is  immediately  apparent  that  the  (real)  eigenvalues  of  the 


linearized  kinetic  matrix  M  are  the  relaxation  times  of  the  coupled 


chemically  reacting  system 


(6.11) 


so  that 

<-  21  “ 

^  1  -t-lW-Tv. 


(6.12) 


We  next  demonstrate  that  Since  M  =  L  ,  one  has 

form  Eq.  (6.1)  after  multiplying  both  sides  of  the  equation  by^'^^: 
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Thennodynainic  stability  requires  that  L  is  symmetric,  i.e. 

and  it  is  an  easy  matter  to  show  from  Eqs.  (2.6)  and  (3.6)  thati/|-is 

symmetric.  Accordingly,  it  follows  that  ^^s  from 

*1/2 

Eq.  (6.3)  after  multiplying  both  sides  of  the  equation  byvil-  ; 

(6.14) 


It  follows  from  Eqs.  (6.13)  and  (6.14)  (the  matrix  ^  has  no  zero  eigen 
values)  that 


(6.15) 


We  may  write 


-I 


JL 


(6.16) 


where  in  the  second 
Si  nue 


rl.T 


relation  we  have  used  the  fact  that  [SlL  ] 


(6.17) 


we  may  substitute  Eq.  (6.16)  to  demonstrate  that 
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a 


VI 


Thus  Eq.  (6.12)  becomes 


I  4  / 


(6.18) 


In  sum,  in  order  to  calculate  <jiui)  for  coupled  chemical  reactions 
one  must  know  U*  and  the  eigenvalues  and  eigenvectors  of  .  While 
this  will  not  always  prove  possible  to  do  analytically,  it  can  easily 
be  accomplished  on  a  computer. 

As  an  example  we  compute  the  exactXxl^)  the  mechanism 
considered  by  Mellen  ^  a2.  displayed  in  Eq.  (5.36).  The  matrix  M  is 
simply 

r  (k.+L) 

h  - 

I  ' 


(  kl-b  It-il 
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The  eigenvalues  are  found  from  the  equation  det  =  det  [JI-m’’’]  -  0; 

(6.20) 


The  tv^o  roots  will  occur  in  the  vicinity  of)c'^lK  +  ^'i')  and  ~  ^  Iz-i). 

A 

If  we  are  interested  in  the  second  root  and  ( kvd- 


the  first  reaction  step  is  fast  compared  with  the  second,  one  is 


justified  in  dropping  the  x  in  the  uppoer  left  hand  corner  of  Eq.  (5.42). 
This  yields  a  single  eigenvalue  which,  as  expected,  is  identical  to 
Eq.  (5.42),  for  the  case  of  a  single,  slow  rate  determining  step 

considered  in  Sec.  V.C. 

The  eigenvalues  are 


A 


Wl,  llAdL 


•4  Jl 

-  X. 


(6.21) 


where 


(6.22) 


The  eigenvalues  are  plainly  real  and  positive. 
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Note  that  for/M-j  >'>4^2’  ^  approaches  the  relaxation  time  Eq.  (5.42). 
The  eigenvalues  of  the  matrix 


r  ^1 

h  --  I 

- 


may  be  expressed  in  terms  of  the  parameters 

IJ-  ^  ,  4 

-  /Vvi, 


r  ^ 


(6.23) 


(6.24) 


so  that 


(6.25) 


and 


cU:  t  -  ^  +  X  =  ^  hi 


(6.26) 
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The  right  hand  eigenvectors  are 


1  ^ 

(6.27a) 

(  i  ) 

corresponding  to  and 


(6.27b) 


corresponding  to  < 

The  left  hand  eigenvectors  are 


(6.28a) 


corresponding  to 


5. 


and 


(6.28b) 


corresponding  to 


Note  that 


(6.29) 
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Normalization  of  these  eigenvectors  is  accomplished  by  requiring 

‘-\'r  ,1,^ 


;  y  ■  ty  •  ( +- )  ■  t-  ■  ^ 


so  that 


and 


cJ-  =  4  i 


According  to  Eq.  (6.10) 


(6.3na) 


(6.30b) 


(6.31a) 


A'% 


and 


1 


(6.31b) 
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Thus  we  have  found  the  complete  solution  to  the  reaction  mechanism 
Eq.  (5.36).  The  frequency  dependent  compressibility  is  [from  Eq.  (6.18)] 
simply 

i  +  iiur, 

where  b+  is  given  by  Eq.  (6.31)  and'J;."^+  =  ^  is  given  by  Eq.  (6.21). 

It  is  evident  that  much  more  complicated  reaction  schemes,  such 
as  those  which  characterize  chemical  relaxation  in  the  ocean,  can 
easily  be  treated  by  the  method  exctly,  if  use  is  made  of  a  computer 
to  determine  the  eigenvectors  and  eigenvalues  for  a  particular 
reaction  mechanism  with  known  rate  coefficients.  Accordingly,  it 
v/ould  be  possible  and  desirable  to  attempt  to  fit  the  entire  spectrum 
of  sea  water  absorption  exactly  in  order  to  verify  the  accuracy  of 
existing  assumed  chemical  reaction  models. 


V 
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VII  Sound  Absorption  from  Systems  Far  From  Equilibrium 

For  systems  near  chemical  equilibrium,  thermodynamics  assures 
that  the  relaxation  times  (eigenvalues)  encountered  are  real  and 
positive.  However,  in  a  steady  state,  far  from  equilibrium  a  coupled 
non-linear  reaction  scheme  may  exhibit  more  complex  behavior. 
Fluctuations  from  the  steady  state  may  be  stable  [Re(A^)>0],  unstable 
[Re(A^)<0],  or  marginally  stable  [Re(x^)=0].  In  addition  fluctuations 
may  exhibit  oscillatory  behavior  [Im(A^)?iO]. 

In  recent  years  there  has  been  considerable  interest  in  the 
behavior  of  such  non-linear  chemical  systems  far  from  equilibrium^ 

In  this  section,  we  wish  to  illustrate  briefly  the  phenomena  that  may 
be  encountered  when  sound  is  absorbed  by  a  chemical  system  in  a 
steady  state.  The  example  discussed  is  not  comprehensive  but  displays, 
in  a  simple  manner,  the  features  that  can  arise  in  the  sound  absorption 
spectrum. 

Accoustic  effects  in  such  systems  have  been  studied  by  Ross 

and  co-workers The  method  of  analysis  followed  here  is  quite 

(Bl 

similar  to  the  method,  discussed  by  Deutch  ^^:,Vequired  to  describe 
the  spectrum  of  light  scattered  from  chemically  reacting  systems  far 
from  equilibrium. 

(91 

A.  The  Prigogine-Lefever  Mechanism^  ' 

Consider  the  model  reaction  scheme  (without  back  reaction) 

k  - =5  X 


(7.1) 
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For  simplicity  in  analysis  we  assume  that  all  the  rate  coefficients 
are  equal  and  that  A=l.  Under  these  circumstances  we  have 


4.  kU, 


where  we  have  introduced  the  two  progress  variables 

T  \jiy  , 


(7.6) 


We  know  that  the  relaxation  times  of  this  system  will  be 
related  to  the  eigenvalues  of  the  matrix  which  takes  the  form 


n 


le 


IS 


(7.7) 


The  two  relaxation  times  which  arise  are 


(7.8) 
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For  B>2,  both  relaxation  times  will  have  a  negative  real  part  since 
Re(Xf)<0.  For  this  case  the  fluctuation  will  grow  exponentially  from 
the  steady  state  signalling  that  the  state  is  absolutely  unstable  and 
hency  physically  unrealizable. 

For  0<B<2,  the  relaxation  time  will  have  positive  real  parts, 
indicating  that  the  candidate  steady  state  is  stable.  However,  in 
this  regieme,  the  relaxation  times  become  complex  with  relaxation 
times  given  by 

(7.9) 


The  imaginary  part  of  the  relaxation  time  indicates  that  the  concen¬ 
tration  variables  perturbed  from  their  steady  state,  for  example  by 
passage  of  a  pressure  disturbance,  will  return  to  equilibrium  in  an 
oscillatory  manner. 


where  '  and  denote  the  real  and  imaginary  parts  of  the  eigen¬ 
values  of  X  respectively. 

It  is  worthwhile  remarking  that  as  the  point  of  marginal 
stability  is  approached  8-^(2-),  the  example  displays  a  slowing  down 
in  the  damping  of  the  chemical  mode  that  one  might  guess  would 
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be  manifest  in  the  spectrum  a(w)  as  a  sharp  peak.  However,  it  must 
be  remembered  that  the  treatment  presented  here  has  ignored  other 
dissipative  transport  effects,  in  particular,  diffusion.  Sufficiently 
near  the  point  of  instability  it  is  necessary  to  take  these  effects 

fq) 

into  account  '  ,  i.e.  one  cannot  rest  the  analysis  on  the  basis  of 

homogeneous  perturbations. 

B .  Significance  of  Complex  Chemical  Relaxation  Times 

The  preceeding  analysis  indicates  that  for  chemical  systems  in 
a  steady  state  far  from  equilibrium,  one  of  the  chemical  modes  will 
possess  a  complex  eigenvalue  in  its  fluctuation  behavior.  According 
to  Eq.  (7.10)  the  equation  of  motion  for  this  mode  is  given  by 

(7.11) 

where^(.  is  equal  to  the  imaginary  part  of  the  eigenvalue  characterizing 
the  mode  and  X  is  equal  to  the  real  part  of  the  eigenvalue. 

It  is  apparent  from  Eq.  (7.11)  that  the  chemical  mode  exhibits 
wave  like  behavior.  This  chemical  wave  is  characterized  by  frequency 
and  damping  Y  .  In  contrast,  when  the  chemical  system  is  at  equili¬ 
brium  and  not  in  a  steady  state4^^=0. 

The  consequences  of  this  type  of  behavior  on  the  frequency 
dependent  compressibility  K:j((a))  can  be  seen  from  Eq.  (4.13). 

In  place  of  the  conventional  expression 
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Xt  ^  "  -^7  ^ 


(7.12) 


one  finds 


^  y.  itJc  ~j 
)r  ^  Y-t  \ . 


(7.13) 


One  finds  that  <y((ja)  contains  frequency  shifts  of  +u)^.  The  chemical 
wave  has  effectively  added  (or  subtracted)  the  frequency  of  the 
chemical  wave  from  the  frequency  of  the  imposed  pressure  wave  w. 

The  pressure  relaxation  in  this  case  may  be  regarded  as  a  sort  of 
Raman  process  with  the  incident  frequency  w,  scattered  into  frequencies 

The  complex  part  of  the  compressibility  <j"((i))  will  be 


f  - - - -  -f  — - - 


so  that  the  Lorentzian  form  for  the  absorption  coefficient  [a(a))/u] 
Eq.  (4.17)  is  replaced  by  the  shifted  Lorentzian 
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oiit^)/)  -  ea 


(7.14) 


in  the  limit  of  low  dispersion.  It  is  dear  from  Eq.  (7.14)  that  the 
absorption  spectrum  is  shifted  by  frequency  as  a  result  of  the 
chemical  waves. 

The  simple  example  considered  here  clearly  indicates  that  sound 
absorption  is  an  important  potential  tool  for  studying  chemical  systems 
in  steady  states  far  from  equilibrium. 
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